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ABSTRACT 

We present a general statistical analysis of Gamma Ray Bursts embedded in 
a host population. If no host generates more than one observed burst, then we 
show that there is a model independent lower bound on the number of hosts, 
H, of the form H > cB^, where B is the number of observed bursts, and c is a 
constant of order one which depends on the confidence level (CL) attached to the 
bound. An analysis by Tegmark et al. (1996) shows that the BATSE 3B catalog 
of 1122 bursts is consistent with no repeaters being present, and assuming 
that this is indeed the case, our result implies a host population with at least 
H = 1.2 X 10^ members. Without the explicit assumption of no repeaters, a 
Bayesian analysis based on the results of Tegmark et al. (1996) can be performed 
which gives the weaker bound of if > 1.7 x 10^ at the 90% CL. In the light of 
the non-detection of identifiable hosts in the small error-boxes associated with 
transient counterparts to GRBs, this result gives a model independent lower 
bound to the number of any rare or exotic hosts. If in fact GRBs are found 
to be associated with a particular sub-class of galaxies, then an analysis along 
the lines presented here can be used to place a lower bound on the fraction of 
galaxies in this sub-class. Another possibility is to treat galaxy clusters (rather 
than individual galaxies) as the host population, provided that the angular size 
of each cluster considered is less than the resolution of the detector. Finally, 
if repeaters are ever detected in a statistically significant manner, this analysis 
can be readily adapted to find upper and lower hmits on H. 
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1. Introduction 

Recently there have been exciting developments in the observational investigation of 
the nature of Gamma Ray Bursts (GRBs) and the hosts in which they may originate. 
(Several excellent reviews of GRBs are available; see e.g. Blaes [1994]; Piran [1995].) For 
GRB970228, an afterglow was detected in the X-Ray and optical wavebands (Van Paradijs 
et al. 1997), and Metzger et al. (1997) have identified absorption hnes in the spectra of the 
optical transient associated with GRB970508. Their observations put the transient beyond 
a redshift of 2; = 0.835. If confirmed, this would establish that GRBs are cosmological. If 
so, then this opens the door to their use as tools for probing cosmological issues such as 
structure formation (Piran & Singh 1997). 

The latest observations herald a new era of more precise, useful information about 
GRBs. The counterparts in the X-ray and optical wavebands allow one to get among other 
things a more precise location of the bursting source. Further, one can perform a detailed 
and intensive survey of the neighboring region of space to seek out the host population 
in which the GRBs are embedded. So far, in spite of exhaustive searches, statistical 
analysis, and the accompanying debate (Fenimore et al. 1993; Larson, McLean & Becklin 
1996; Schaefer et al. 1997; Blaes et al. 1997; Kolatt & Piran 1996; Hurley et al. 1997), 
GRBs have not clearly and unambiguously been associated with any other well-known 
population of astronomical objects. In particular, no known or identifiable galaxy has been 
associated with either of the two GRBs for which transient counterparts have enabled a 
small error-box. Regarding repetitions of GRBs from the same host, while there was an 
early claim of evidence for this by Quashnock & Lamb (1993), their analysis was called 
into question by Narayan & Piran (1993). More recently, Tegmark et al. (1996) concluded 
that the data is consistent with no repetition of GRBs. Thus all of the recent data has not 
yet resolved the mystery of the nature of the host population in which GRB sources are 
embedded. 

It is therefore important to deduce all possible constraints which can be placed on 
the host population; model-independent constraints are particularly valuable. Obviously, 
there must be at least as many hosts as there have been GRBs, assuming no repeaters. 
This is clearly not a significant constraint; the total number of galaxies is about 10^, 
corresponding to about 10~^ bursts per galaxy per year. It is possible, however, that the 
host population consists of exotic or otherwise rare cosmological objects. In this Letter, 
we perform a statistical analysis which demonstrates that, if a total of B bursts have been 
seen, and we assume that there are no repeaters, then the number of hosts H is likely to 
exceed cB'^, where c is a numerical constant which depends on the desired confidence level 
and the precise form of the statistic (e.g., choice of the prior distribution in a Bayesian 
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analysis) , but not on any details of the GRB model. If we do not make the assumption of 
no repeaters, but instead appropriately weight the possibilities based on observational data, 
then the lower bound is reduced to i^min ~ B"^ / {^2)1 where (A^2) is the expected number of 
repeaters in the data. In this analysis, a "host" is any object whose angular size is much 
less than the resolution of the detector. 

The BATSE data is consistent with no repeaters in the analyzed catalog of 1122 bursts. 
An analysis by Tcgmark ct al. (1996) shows that the number of repeaters N2 is less than 10 
at the 95% CL. If wc assume that there are in fact no repeaters in this catalog, then we can 
conclude that the the population of hosts from which the GRBs may be originating should 
contain at least H = 1.2 x 10^ members. This implies that about 0.1% of the total number 
of galaxies could be hosts, if GRBs are embedded in galaxies. This could be a significant 
constraint if GRBs are found to be associated with a specific sub-class of galaxies. This 
constraint would also be significant if ongoing searches for potential objects within the error 
boxes of the transients fails to turn up identifiable normal galaxies or other well-known and 
abundant astronomical objects. Furthermore, if repeaters are ever detected in a statistically 
significant manner, our analysis can be readily adapted to find upper and lower limits on 
the number of hosts H. 

We now turn to the statistical analysis of this problem which allows us to draw the 
advertised conclusions. 



We suppose that there is a population of H hosts, each of which is producing burst 
bursts at an individual rate of 7^, i — 1, . . . , if, and that each of these hosts has been 
observed for an amount of time tj. Then the number of bursts which have been observed 
from each host is expected to be Ai = ^iti. Call the actual number of observed bursts from 
each host B^, and the total number of observed bursts B = J2iLi A- 

We assume that the bursts from each host are uncorrelated, so that the probability of 
observing B^ bursts from host i is given by a Poisson distribution 



2. The General Statistical Analysis 
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Then the probability that a total of B bursts are observed is 
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(AH)^ exp{-AH) 

B\ ' ^ ' 

where A is the average value A = {l/H) E,=i Ai. 

We wish to calculate that probability that B bursts are seen, but that there are no 
repetitions; this is given by 

P(5,norep|{A,},//)= ^ ... ^ 5b,+...+Bu,b P{BMi) ■ ■ ■ P{Bh\Ah) (3) 

Bi=0 B„=0 

To simplify the calculation, we first treat the idealized situation in which Ai = A for each 
host. In this case, eq. (0) becomes 

P{B, no repi A, H) = A^ exp(-A/7) , (4) 

The probability that no repetitions are seen, given that B bursts are observed, is given by 
the ratio 

F [no Tep\B, A, H) = ■ 



PiB\A,H) 
H\ 

(5) 



H^{H-B)\ ' 

which we see is independent of A. Assuming H ^ B ^ 1, and using Stirling's formula for 
the factorials, we get 

P{noTep\B,A,H) = exp{-B^/2H) . (6) 

This is our key result. It tells us that we are unlikely to have seen no repetitions among 
B bursts if H is significantly smaller that B^. This can be made more precise by using a 
Bayesian analysis; see, e.g., Loredo and Wasserman (1995). Assuming that no repetitions 
are seen, the probability that the number of hosts is between H and H + dH is proportional 
to the likelihood P(no Tep\B, A, H) times a prior distribution; if we make the standard 
scale-invariant choice dH/H with an upper cutoff at i^max = 10^, we find that H > B^ with 
a confidence level of approximately 90% (the exact value depends weakly on B and -ffmax)- 

If some hosts generate repeat bursts, the analysis above is changed somewhat. Let N^y 
be the number of hosts which generated u bursts each; then H = Nq + Ni + N2 + . . . and 
B = N1 + 2N2 + 3N3 + .... We assume A^o > A^i > ^^2,3,...- Then eq. (g) is replaced with 

HI B\ 

= ^ ; exp{-By2H) exp{-2Ni/B)6N,oSN,o • • • • (7) 
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Terms with A^^, 7^ for u > 3 are suppressed by at least one factor of B/H. Henceforth, we 
set = N4 = . . . = 0. 

Let us now consider what can be deduced from our formulas when they are applied to 
the BATSE 3B catalog oi B = 1122 bursts. First we will examine the results of Tegmark et 
al (1996) concerning repeaters. To gain an analytic understanding of these results, we will 
treat the simplified case in which the beam- width a is the same for every burst; then their 
statistic R becomes 

^=^EEexpMj/4a2), (8) 
^-'^ i=i i=i 

where 6ij is the angle on the sky between bursts i and j. 

We now compute the expected values of R and R^, assuming that the hosts are 
uniformly distributed on the sky, that the angular size of each host is much less than a, and 
that each the observed location of each burst is displaced from its true value by a random 
angle selected from a gaussian distribution of width a. (The latter two points are relevant 
only for bursts from repeaters.) In the flat-sky approximation (a ^ 1), we find 

{R) = \Ba' + \f (9) 
(Ai?)' = {R^)-{RY 

= T-y + ieB-'f, (10) 

where / = 2N2/B is the fraction of bursts which came from repeaters. If each host has an 
angular size 0, then the coefficient of / in eq. (^) is reduced; if the burst probability across 
a host has a gaussian profile, then the reduction factor is a^/ (cr^ + 0^). Not including such a 
factor in our subsequent analysis implies that each of our hypothetical hosts has an angular 
size which is much less than a. 

Eqs. (p|JTO|) can be understood heuristically. Roughly speaking, {R) is the expected 
number of bursts which lie within an angular distance a of any one particular burst. There 
is a contribution to (R) due to bursts from different hosts; this contribution is of order 
Ba^. There is a second contribution if the particular burst in question is a repeater; the 
probability that this is the case is of order /. This explains the structure of eq. (P). To 
understand the variance (Ai?)^, we note that the Ba^ term in (R) comes from a Poisson 
process, so its variance is also Ba"^. However, averaging over the particular burst in question 
reduces the variance by a factor of B. The term of order / in (Ai?)^, which in practice is 
not important, can be similarly understood. 

The probability distribution for R should be approximately gaussian, 

PiR\N2,B) = ^^exp[-(i?- (i?))V2(Ai?)2] , (11) 
V 27rArt 
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where i?av and AR are given in terms of A^2 = fB/2 and B by eqs. (^,0). For BATSE, in 
which each burst has its own value of a, and in which sky coverage is not uniform, eqs. (p|,p!OD 
should be approximately correct, but possibly with different numerical coefficients. For 
this case, we turn to the Monte Carlo results of Tegmark et al. (1996). Their Figure 1 
shows P{r)dr vs. R for several different values of /. From this, we can read off 
(R) = 0.316 + 0.15 / and AR = 0.0068 for / = 0. This is consistent with our eqs. (g^) 
with a value of a ~ 2°. (We neglect the / dependence of AR, which is difficult to extract 
accurately, since it will not be important to our conclusions.) The experimental value of 
R for BATSE is -Rexp = 0.308. Assuming that there are in fact no repeaters, then the 
probability that the BATSE data would yield a value of R which is less than -Rexp is 0.11. 
Eq. ( pUD is also consistent with the results of Tegmark et al. (1996) that / < 0.018 at 95% 
CL and / < 0.049 at the 99% CL. 

We now combine eq. (0) and eq. (pAf) to get 

B/2 

P(i?exp|5, A,H)= J2 PiRe.p\N2, B)P{N2\B, A, H) . (12) 

N2=0 

Again performing a Bayesian analysis for H with the scale-invariant prior distribution 
dH/H and an upper cutoff of i/max = 10^, we find H > 1.7 x 10^ with a CL of 91%; this is 
smaller than B^ by a factor of 7.4. This factor is roughly the expected number of repeaters 
in a Bayesian analysis, 

^""'^ - E^,.o^(^exp|Ar„5) ' ^^2) 



where P{Rexp\N2, B) is given by eq. ([Tl|), and we have implicitly assumed a uniform prior 
distribution for A^2; in the present case, eq. (0) yields {N2) = 12. 

Turning to future experiments, we first recall that BATSE happened to produce a low 
value of -Rexp- A larger catalog for BATSE which instead yielded the most likely value when 
/ = 0, i?exp = {R)f=o = \Ba\ would result in P(i?exp|/) « exp[-(0.15 /)V2(AP)2] for 
/ > 0. Assuming (for example) a catalog with B = 3000 bursts, we find that H > 3.7 x 10^ 
at the 90% CL. More significant improvements in the lower bound on H would require an 
instrument with a better angular resolution. 

We now turn to the more general case of a population of hosts, each with an individual 
Ai (the burst production rate times observation time for that host). We will make the mild 
assumption that the distribution of Aj's has moments 

^J:A^^c^A- (14) 
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which remain finite in the limit of large H. Note that Ci = 1 and C2 > 1. From here on we 
will concentrate on the simplest case of no repetitions. 

The probability that we see B bursts with no repetitions is given by eq. (||). We can 
write the Kronecker symbol as 

Sb.^...^b,,b z-^z^^ . . . z^- , (15) 

where the contour encloses the origin. Each sum in eq. (j^) now yields a factor of 
exp(— 74j)(l + zAi), and so we find 

/dz ^ 
zmz 

/dz \ ^ 1 
exp Vlog(l + zA) -51og^- A/J . (16) 
2'Kiz \~i J 
1=1 



We now treat both H and B as large, and evaluate the integral by steepest descent. Let 
F{z) be the argument of the exponential in the second hne of eq. (^); the point Zq of 
steepest descent is then given by F\zq) = 0, or 

Ai B 

= — —A 

j=l 1 + ZoAi Zq 

= H[A-zoic2A^) + z^,ic,A^) + ...]-Bz,K (17) 
Taking B <^ H, this can be solved by power series to yield 

ZoA = (B/H) + C2{B/Hf + {2cl - c^){B/Hf + .... (18) 
The probability of getting B bursts with no repetitions is now given by 

P(5,no rep|{A,},/7) = [27rF"(^o)]"'/'^o"' exp F(^o) • (19) 

After dividing by P{B\A, H) we find that the probability that no repetitions are seen, given 
that B bursts are observed, is 

P(no rep|5, {A,}, H) = exp{-C2By2H) (20) 

in the limit H ^ B ^ 1. This is qualitatively the same as our previous result, eq. (|]), 
except that now the relative value of the second moment of the distribution enters as well. 
The coefficient C2 must be greater than one, so this correction can only increase the lower 
bound on H. 
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3. Implications for the hosts of GRBs and Discussion 

We now turn to the imphcations of the preceding analysis to the specific case of GRBs 
and the host population in which the GRBs may be embedded. We have seen that if we 
observe B GRBs, which originated from a host population of H hosts, and further if assume 
that no host generated more than one burst, then we can conclude that H must be bounded 
below by 5^, at a confidence level of roughly 90%. If a probability- weighted average over 
the number of repeaters is taken instead, then the lower bound on H drops to a number 
which is of order B"^/ {N2), where {N2) is the expected number of repeaters. This is given 
by eq. (|13|), which would typically yield (A^2) ~ Bcr (and hence -f^min ~ B/cr), where a is 
the angular resolution of the detector. Our constraint on H can be significant for models in 
which the host population for GRBs consists of exotic objects of any kind, such as galaxies 
with some particular morphological feature or other specific signature such as a high star 
formation rate. Another possibility is to consider clusters of galaxies as hosts (rather than 
individual galaxies) , provided that the angular size of the clusters considered is much less 
than a. 

If repeaters are ever detected statistically through a value of -Rexp which is larger than 
(R) by several AR, this analysis can be readily adapted to find upper and lower limits on 
the number of hosts H. The most likely value will be if = B^/ {N2), where {N2) is the most 
likely value of the number of repeaters. 

Finally, we restate our main point: simply by knowing that there are ~ 10^ bursts with 
no sign of repeaters, we can immediately infer a model-independent lower bound on the 
number of hosts in which GRBs may be embedded which is larger by more than two orders 
of magnitude. 

We thank Omer Blaes for many helpful discussions, and the referee for several 
important comments which improved the paper. This work was supported in part by NSF 
Grant PHY-91-16964. 
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